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ABSTRACT 

We obtained self-similar solutions of relativistically expanding magnetic loops taking 
into account the azimuthal magnetic fields. We neglect stellar rotation and assume ax- 
isymmetry and a purely radial flow. As the magnetic loops expand, the initial dipole 
magnetic field is stretched into the radial direction. When the expansion speed ap- 
proaches the light speed, the displacement current reduces the toroidal current and 
modifies the distribution of the plasma lifted up from the central star. Since these self- 
similar solutions describe the free expansion of the magnetic loops, i.e., Dv/Dt = 0, 
the equations of motion are similar to those of the static relativistic magnetohydro- 
dynamics. This allows us to estimate the total energy stored in the magnetic loops 
by applying the virial theorem. This energy is comparable to that of the giant flares 
observed in magnetars. 
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1 INTRODUCTION 

Soft gamma-ray repeat ers (SGRs) are believed to be a young neut ron star with strong magnetic fields (~ 10^^ G), namely 
the magnetar (see, e.g.. Woods &: Thompson 20061 : Mereghettil 20081 . for review). The magnetic fields inside the magnetar are 
amplified by the dynamo mechanism at the birth of the neutron star. The Lorentz force stressing the crust of the magnetar 
balances with the rigidity of the crust. When the critical twi st is accumulated, the magn etic twist injected into the magnetar 
magnetosphere will trigger the expansion of magnetic loops ([Thompson fc Duncan 2001 ). The magnetic re connection taking 
place inside the expanding magnetic loops can be responsible for SGR flares ( Woods et al. 2001 ; Lvutikov| [2006) . 

Recently, relat i vistic simulations have been performed to study the dy n amics of th e magn etospheres of neutron stars 



( Komissarov 



relat i vistic simulations nave Deen periormea to stuay tne gy n amics or tn e magn etospneres ot neutron star; 
2OO2I : lAsano et all I2OO5I : ISpitkovskvl bood : iKomissaro-JI bood ). ISpitkovskvl (|2005h reported the results of 2 



dimensional relativistic force-free simulations of the magnetar flares triggered by the injection of the mag netic tw i sts at 
the footpoints of the loops. When the critical twist is accumulated, the magnetic loops expand relativistically. Asanol (2007) 
carried out 2-dimensional relativistic force-free simulations of ex panding magn etic loops and showed that the Lorentz factor 
deflned by the drift velocity Vd ~ c{E x B)/B'^ exceeds 10 (see. IUchidalll997l . for the deflnition of the drift velocity). These 
simulations indicate that the magnetic loops e xpand self-similarly. 

Assuming relativistic force-free dynamics , Lvutikov fc BlandfordI ( 20031 ) obtained self-similar solutions of the spherically 
expanding magnetic shell. iPrendergasti (|2005l ) found self-similar solutions of the relativistic force-free field. In these studies of 
force-fre e dynamics, gas pressure and inertial terms are neglected. In the framework of the relativis tic m a gneto hydrodynamics 
(MHD) ■ iLvutikovl (|2002l ) found self-similar solutions of the spherically expanding magnetic shells. iLowl (|l982l ) obtained non- 
relativistic self-similar MHD solu t ions o f the expanding magnetic loops in solar fiares or supernovae explosion by assuming 
axisymmetry. Subsequently, iLowl (|l984l ) extended his model to the case in cluding toroidal m agnetic fields and applied it 
to solar coronal mass ejections (CMEs). The latter model was employed bv lStone et al.l (|1992l ) as a test problem to check 
the validity and accuracy of ajcisymmetric MHD codes. In magnetar flares, the magnetic loops may be twisted by the shear 
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motion at the footpoints of the loops. The shear motion generates Alfven waves propagating along the field lines. Such twisted 
magnetic loops expand by the enhanced magnetic pressure by the toroidal magnetic fields. Thus we should include the toroidal 
magnetic field to study the evolution of magnetic loops during magnetar flares. Also the relativistic effects should be included. 
The characteristic wave speed in the magnetar magnetosphere approaches the light speed because of the strong magnetic 
fields. Thus our aim is to obtain relativistic self-similar MHD solutions of e xpanding m agnetic loops taking into account the 
toroidal magnetic fields by extending the non-relativistic solutions found by Low ( 19821 ) . 

This paper is organized as follows; in § [2l we present the relativistic ideal MHD equations and introduce a self-similar 
parameter which depends on both radial distance from the centre of the star and time. In § |3l we obtain self-similar solutions. 
The physical properties of these solutions are discussed in § [J] We summarize the results in § [S] 



2 SELF-SIMILAR MHD EQUATIONS 

In the following, we take the light speed as unity. The complete set of relativistic ideal MHD equations is 

^(7P) + V- (7P1^) =0, (1) 

P7 + ■ V)] (?7«) = -Vp + p.E +jxB- ^^^er, (2) 



lniL)=0, (3) 



\/-E = 4npe, (4) 
V ■ B = 0, (5) 

^+VxiS = 0, (6) 

^=VxB-4^j, (7) 
E = -V X B, (8) 

where E , B , j ,v,j, pe, p,p,T are the electric field, the magnetic field, the current density, the velocity, the Lorentz factor, the 
charge density, the mass density, the pressure and the specific heat ratio, respectively. The vector Br is a unit vector in the 
radial direction. We include the gravity by a point mass M as an external force. Here G is the gravitational constant, and r 
is the distance from the centre of the star. The relativistic specific enthalpy ^ is defined as 

where e is the energy density of matter including the photon energy coupled with the plasma. In SGR outbursts, since the 
luminosity much exceeds the Eddington luminosity, radiation energy density can exceed the thermal energy of the plasma. In 
the following pressure p includes the contribution from the radiation pressure. 

In this paper, we consider relativistic self-similar expansions of magnetic loops which started expansion a,t t — ts by loss 
of dynamical equilibrium and entered into a self-similar stage at t = to > is. We do not consider the evolution of the loops 
before t = to- 

For simplicity, we ignore the stellar rotation and assume axisymmetry. We can express the axisymmetric magnetic field 
in terms of two scalar functions A and B as 

(10) 



r dd ' dr 

in the polar coordinates {r,9,(f>). The scalar function A{t,r,9) denotes the magnetic fiux, whose contours coincide with 
magnetic field lines projected on to the r — 6 plane. 

We further assume that the fiuid flow is purely radial; 

V ^ v{t,r,e)er. (11) 

Equations ([T}, ([2]), ([Sjl, and (O are then expressed as 

^ + 4^fe^=0, (12) 
at r-^ or 
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A 2 ■ 2 adp ,^ 2nd 

47rr sm 0— + (1 — t; )B 



de 



^ ' dr de de 



— +v — 

dt dr 



In- 



dB dA 

'de^'de 

dB ^ d{vB) 
dr dt 

= 0, 



d f dA 



dt V dr 



dAdv 

"S— — = 0, 
dr do 



dA ^dA _ ^ 
dt dr ' 
dB ^ djvB) _ 
dt dr ' 

where we used the MHD condition given by ^ and introduced the operator 



d^ 



(r.e) 



+ ■ 



sine d 



( — -) 



de 



(14) 

(15) 

(16) 

(17) 
(18) 

(19) 



Since our aim is to obtain self-similar solutions of these relativistic MHD equations, we assume that the time evolution is 
governed by the self-similar variable: 
r 



(20) 



where Z{t) is an arbitrary function of time. We further assume that the flux function A depends on time t and the radial 
distance r through the self-similar variable rj, as 



Ait,r,e) = A{y,e). 

When equation (|21|l is satisfied, the radial velocity v has a form 
i; = r]Z, 



(21) 



(22) 



from equation (|17|l . Here dot denotes the time derivative. Equation (|22l) implies that the radial velocity v does not depend 
on the polar angle e. It then follows from equations p2p and psp that 



p(t,r,e)^it,r) = Z-^t)D{r^,e), 
B{t,r,e) = Z-\t)Q{r,,e), 



(23) 
(24) 



where Q and D are arbitrary functions of ri and e. These relations indicate that the magnetic flux and the total mass are 
conserved. Next we take the pressure p as p(t,r,e) = Z''P{rj,e). Substituting this equation into equations (|14p and ((TSJ, we 
obtain 



. 2rjl+4, ■ 2 ^dP dA 



C,,M+{VZZ^2,Z^)'J-~,^Z^'^ 



+ (l-^^Z^)O^=0, 



Tri^ZZ 



+ (3r-f0 = 0, 



1 - 772^2 

where we introduced an operator L^^i.e)' 



d^ 



svae d 



V sine 961 ^2'^('-.'')- 



(-7.9) - Q^2 ^ ^2 QQ 

To satisfy these equations, p, Z and F should have forms 

p{t,r,e)^z'^p{n,e), 
Z{t) = t, 
and 

This adiabatic index corresponds to the radiation pressure dominant plasma. Thus our model can describe the evolution of a 
fireball confined by magnetic fields. 

Equations (|23|l and p8|) indicate that the magnetic loops expand adiabatically. By using equations (|22|) . (|23|) . (|24|) . (|28p . 
H29[l and (|30|) . equations p3|) . (|15p and (|25[) are expressed as 



(25) 
(26) 

(27) 

(28) 
(29) 

(30) 
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Figure 1. The flux function A as a. function of r) for the dipolar solution (solid curve), the shell solution (dashed curve), and the flux 
rope solution (dot-dashed curve). The parameters a and b denote the outer and inner boundaries of the shell, respectively. 



dP 1 
drj Atvi]^ sin^ i 



dA 



2.dAdQ OA d r,, 2.^1 „ 



, 2 . 2^dP ^dA 



^{ri,e)A - 



d 



drj \^ drj 

From equation (|32p . a formal solution of Q is obtained as 



,dA 



(31) 
(32) 
(33) 

(34) 



1-772' 

where Q is an arbitrary function. 

Self-similar solutions can be constructed as follows. First we prescribe an arbitrary function A{rj,6) (or Q{rj,6)). Then, 
equation (I32|) determines the function Q{rj,9) (or A(rj,6)). Functions A and Q determine the pressure P{rj,9) according to 
equation (I33|) . Finally, the density function D{rj,9) is obtained by equation (|3ip . 

Note that from the equation (|22p and (|29p . the radial velocity has a simple form as 



r 

V = -. 

t 



(35) 



Since the time derivative of the velocity becomes zero, i.e., Dv/Dt — 0, equations (|3ip - p3p describe the freely expanding 
solution. This means that there is a reference frame that all forces balance. By substituting equations (I23p . (|24p . (|28|l . (|29p . 
and p5|) into the equations of motion ([2]), we obtain 



r j^v'^p 



Vp + peB + J X B = 0. 



(36) 



r - 1 r - - 

The first term on the left hand side comes from the inertia. For convenience, we call this term as a thermal inertial term 
throughout this paper. When we neglect the terms of order (w/c)^, equation (|36[> reduces to the equations of the force balance 
in non-relativistic MHD. 



3 SELF-SIMILAR SOLUTIONS 

In the previous section, we derived relativistic self-similar MHD equations, (|3ip . (I32p and (|33p . In this section, we obtain 
solutions of these equations by imposing appropriate boundary conditions. As mentioned in the previous section, the toroidal 
magnetic field, the pressure, and the gas density are calculated by assigning the flux function A{rj,6). In the following, we 
introduce three kinds of flux functions and obtain explicit forms of other variables. 



3.1 Construction of Solutions 



We assume that the expanding magnetic loops have a spherical outer boundary at r = R{t). 

A sirnple s olution of the expanding magnetic loops is that the poloidal magnetic field is dipolar near the surface of the 
star (|LowII 19821 ). The magnetic field should be tangential to the spherical surface r = R{t) at all time. Such a solution can be 
constructed by 
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Tj/a rj/a rffa 

Figure 2. Contour plots of the magnetic flux A wtiicti constructs tlie dipolar, shell, and flux rope solutions from left to right, respectively. 
The parameters are taken as a = 0.8, b/a = 0.7, and k = 7r/[4(a — b)]. 



2 2 



where Ao and a are constants. The radius R{t) where ^ = is given by 
R{t) = at. 



(37) 



(38) 



We hereafter call the solution constructed from equation (|37l) as dipolar solution. 

Solid curve in Fig. [1] shows the flux function ^4 as a function of for the dipolar solution. Contour plots of A for dipolar 
solution is shown in the left panel of Fig. [2] 

When the flux function is given by equation p7|l . the magnetic flux crossing the annulus at the equatorial plane 6 — 7r/2 
decreases with radius (see Fig. [1]). In actual MHD explosion, the magnetic flux can be swept up into a thin shell just behind 
the loop top. The shell boundaries are assumed to be at r = bt and r — at (region II, see Fig. [!}. Such a self-similar field can 
be constructed by 

V S b), 



where 
A{v) = 1 - 
T{v) = Hv 



Aqc^ sin^ 9, 
Aoa^Airi) sai' 6, 



(region I : 
(region II 



b <ri^ a), 



sin* T(a) 
-b), 



(39) 

(40) 
(41) 



and a, b and k are constants ( Low IQSj ). The flux functions in region I (77 ^ fo) and region II (b < 77 ^ a) are connected 
smoothly ai rj — b. The loop boundary locates at r = at, where A — 0. 

The fiux function for this solution is shown by a dashed curve in Fig. [T] It can be easily shown that the magnetic field 
lines projected on to the r — 6 plane are all radial in region I. We call the solution constructed from equation (|39p as shell 
solution. The middle panel of Fig. [2] shows the contours of A for the shell solution. 

Another solution is that we call flux rope solution. As the magnetic loops expand, a current sheet is formed inside the 
magne tic loops. It is suggested that the m agnetic reconnection taking place in the current sheet is responsible for the SGR 



flares (Woods et al. 200 



jjg ested that the m agnetic reconnection taJiing place m the current sheet is responsible tor the sLxK 
J Lvutikov 20061 ). When the magnetic reconnection takes place, flux ropes (namely plasmoids) are 



formed inside magnetic loops. The flux function should then have a local maximum inside the fiux rope. Such a solution can 
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be constructed by 

Aoa^ 



n2 



(region I : 77 ^ 6), 



MV,0) = { \/ (42) 



A(77) sin 6, (region II : b < i] ^ a), 



where Aq, a are constants and A(77) is given by equation (|40|l . This function is shown by a dot-dashed curve in Fig.[T] It has a 
local maximum in the domain b < r) < a (see Fig. [T]). The contours of A for the flux rope solution is shown in the right panel 
of Fig. [5] Flux ropes appear behind the shell. 

3.2 Dipolar Solutions 

Dipolar solutions are constructed by the flux function specified by equation (|37|) . The azimuthal magnetic fields can be 
obtained by substituting equation p7|l into equation p2p as 

Qin, e) = Y. If ~fj+^ ^- (43) 

(I — ry^) +4 

where Qo,n are constants. Note that the solutions H37[) and l|43p satisfy the formal solution given by equation (|34|) . Substituting 
equations l|37[) and (|43l) into equation (|33p . we obtain the pressure function P: 

P{V,9) = Poiv) + Pa(v,0) + PQiv,e), (44) 
where -Po(?7) is an arbitrary function arisen from the integration and Pa and Pq are given by 

PAiv, 0) = T%7T^^t (2«' - 3a^r?^ - 77^ + 2,,«) sin^ (45) 




nQo,mQo.n {a^ — rj^) ^ . m+71-2 ^ <■ , 

■ sm 6^, tor m + n 7^ 2, 



47r(m + n — 2) „2|'i _ «2^i+^t^ 
J'Q(^,e) = <! '"+"^^ ^ ^ 2^2^^ (46) 

nQo„nOo,n a -r, .^^^^^^^^^^ for m + n = 2. 



™+„^2 4" r,2(l-^^)^ 
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Substituting equations (|37|l . (|43|l . (|44p . (|45p . and (|46p into (|3ip . the density function D can be determined as 

D{'n,e)^Do{ri) + DA{r,.e) + DQ{'n,0), 



where 



4r?Po d^o 



GM V 1 - J?^ rf?? 



AnGMrj^ (1 - 772)3 

nQo,mQo 



a^(8 - 127?^ + Sry'') - r]^{5 - 677^)] sin^ 9, 



Dq{v,0) = < 



E 

m + Ti 7^ 2 



47rGM(m + n - 2) 



7?(1- 772)2+^ 



Stt 



nQo,mOo,n -J?^(2 - - r;2) + 2(2a2 - a^r;^ - r?*) log(sin6l) 



GM 



77(1-772)! 



for 771 + 77 7^ 2, 



for 771 + 71 = 2. 



(47) 

(48) 
(49) 

(50) 



The parameters 777 and n correspond to the Fourier modes in the 6 direction. These parameters should be determined by 
the boundary condition on the surface of the central star where magnetic twist is injected. 

Equation (|44|) and 1)47^ indicate that the solution consists of three parts, Po, Pa and Pq (or Do, Da and Dq). The 
arbitrary function Poiv) describes an isotropic pressure in the region r < R{t). The isotro pic density profile Do{rj) is related 
to Po through equation (|48p . This equation is similar to that in non-relativistic model ( Low! 1 19821 ). In the non-relativistic 



model, gravity is supported by the gradient of Po(77). In the relativistic case, relativistic correction of the plasma inertia cannot 
be ignored. This effect is included in the first term in the right hand side of equation (|48[) . Other functions Pa and Pq (or 
Da and Dq) come from the interaction with the electromagnetic force. Note that the plasma pressure Pq, which balances 
with the electromagnetic force produced by the toroidal magnetic field, is always negative. This suggests that the pressure is 
smaller for larger toroidal magnetic fields. 

Fig. |3] shows the contour plots of the magnetic fiux A (left), the poloidal part of the pressure Pa (centre), and that of 
the gas density Da (right), while Fig. |4] shows the contour plots of the toroidal magnetic field (left), the toroidal part of 
the pressure Pq (centre), and that of the gas density Dq (right) in the 77/a — 6 plane for 771 = 77 = 4 and a — 0.8. 

The magnetic field is explicitly expressed as 
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,cosge,. + — smge, + )^— sm 0e,, (51) 



where e^, eg, and are unit vectors in r, 9, and directions in the polar coordinate, respectively. Note that Br and are 
zero at r = R(t) but Be is not zero and it depends on time when a ^ 1. We will discuss the physical meaning of this result 
later in i]3.3l 

In later stage, the magnetic field becomes stationary, 
lim B — cos 9er, (52) 

and the magnetic field becomes radial. In the limit t ^ r, the pressure and the gas density inside the magnetic loop are given 
by 

limp=^sin^e+l(Por,*)| (53) 
t—oo nGMr^ ^ / 1.7=0 

Since the toroidal magnetic field tends to be zero in this limit, the pressure and density do not depend on the amplitude of 
the toroidal magnetic fields. 



3.3 Shell Solutions 

Shell solutions are constructed from the fiux function (|39|l . By substituting equation (|39|l into equation (|32p . the function Q 
can be obtained as 



'{r,,e)^Q;,^^, (55) 



fid) 




-l (^\ I I I I l_l I I I I I I I I I l_l I L 

■a.a a:2 (3 a o.fi i.ii 

r}Ut 



4flr(m + 11-2) 



Relativistic Self-similar Expansion 9 
2irGM(m + n - 2} . 



\M\ — I — I — I — I — n — I — pn — I — I — I — I — I — I — [— 1 — I — r 




.\ ^1 I I I I I I I I I I I I I I I I I I I 

0.0 0,2. 0.4 O.t) 0.S 1.0 



TjUi 



-0.5- 



Or 



I ' ' ' I ' ' ' I ' ' ' I 




.[ ()l I I I I I I l—l I I I l_l I I I I I L 

0.0 0:2 0.4 0.6 0.« L.O 



Figure 6. Contour plots of the toroidal magnetic field iJ^ (left), the toroidal part of the pressure Pq (centre), and that of the gas density 
Dq (right) for the shell solution in rj/a — 6 plane when a = 0.8, b = 0.75a, and n = m = 4. 



where f(6) is an arbitrary function of 6, and Qq and Qo^„ are constants. The subscripts I and // mean that the function is 
defined in region I and in region II, respectively. The arbitrary function f{9) can be determined by applying the boundary 
condition that magnetic field should be connected smoothly at 77 = 6, 

Q'iv = b,e) = Q''{7^ = b,9). 



By using the boundary condition, the function / is given by 
/(e) = ^sin"0. 

n 

and the function is obtained as 

n 

The constants Qo_„ and Qo^„ should be related by 
Qo,n = Qo.n = Qo'n sin^" T{a). 



(57) 



(58) 



(59) 



(60) 



from the boundary condition (|57|) . Substituting equations (|39|) . (|56l) . (|59p and (|60p into equation (|33|) . we obtain the pressure 
function P{ri,8). The density function D{ri, 9) is obtained from equation (I31|l . The functions Q, P, and D obtained in region 
I and region II are given in appendix [K\ 

The pressure and the gas density consist of three parts, the isotropic part Po and parts representing the interaction with 
the electromagnetic force by the poloidal and toroidal components of the magnetic field, Pa and Pq , similarly to the dipolar 
solutions (see equations (IA2|l and (|A9[l . for the pressure and equations (|A5P and (|A12[l for the gas density). 

Fig. [5] shows the contour plots of the magnetic fiux A (left), the poloidal part of the pressure Pa (centre), and that of 
the gas density Da (right) in rj/a — 9 plane. Fig. [6] shows contour plots of the toroidal magnetic field (left), the toroidal 
part of the pressure Pq (centre), and that of the gas density Dq (right) in rj/a — 9 plane. The parameters are taken to be 
a — 0.8, b — 0.75a, and m — n — 4 in both figures. A shell structure appears behind the loop top. 

The pressure Pq is always negative (see equations (|A4|l and (|Alip and the middle panel of Fig. |6]) and its amplitudes 
is proportional to that of the toroidal magnetic fields, Qo,n- This indicates that the pressure is smaller for a larger toroidal 
magnetic field. 

The magnetic field is explicitly given by 
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■ cos 9er + Qo,n , „ ^ — 5Tsin" ^ Oe^,, {r]^b), 



r 



2 



^ ^ 2Aoa\ , , AAoa^ , sin^T(r/t)cosT(r/t) . ^ ^ , a , ,,tsin"-^6l „ , , ^^"^^ 

— ^A rA cos6»e, + — ^^fc U^^-^-±J. smOeg +} Qo.nA^ r/f — -e<^, b < 77 S a ■ 

rt sm r(a) r(t^ — r^j 

n 

Similarly to the dipolar solution, the shell solutions have the parameter m which corresponds to the Fourier modes in 
the polar angle 9. These modes and the corresponding amplitude Qo,m of the toroidal magnetic fields should be determined 
by the boundary condition at the surface of the central star where the magnetic twist is injected. 

In contrast to the dipolar solution, the magnetic field lines do not cross the equatorial plane in region I (see the left 
panel in Fig. [5}. Note that in the limit that r, the magnetic fields and plasma distribution approach those of the dipolar 
solution, given by ^E^, ((53]) and (f54)l . 

At the boundary r = R{t), the field components Br and are exactly zero, but Bo is not zero. Since the Poynting fiux 
S — {E X B)/{4:n) is not zero at r = R{t), the energy fiux will be transmitted to the region outside the boundary at r = R{t). 
When cosr(a) — 0, since the magnetic field vanishes at r = R{t), the energy is not transferred to r > R{t). This happens 
when the constant k is given by 

where I is an integer number. 

Fig. [7] shows the distributions of the Poynting fiux ASr = r*7r5r/(Aoa^ sin^ 9) (left panel) and ADa = DAivGMif{l - 
if) / [2 A^a^ siv? 6\ (right panel) for a = 0.8, h/a — 0.75, and B^ — for shell solutions. Solid curve denotes that for k — 
7r/[4(a — b)], while dashed and dot-dashed ones do for k — 7r/[2(a — b)] and k = 37r/[4(a — b)], respectively. When k satisfies 
equation (|62[) . Be(t,r = at, 6) — and thus Sr{t,r = at, 9) — 0. Electromagnetic energy is not transmitted ahead of the loop 
top. When k does not satisfy equation (|62|) . the Poynting fiux Sr at r = at is not zero and the electromagnetic energy is 
transmitted to r > at. The physical interpretation of the condition given in equation (|62|) is as follows. 

Let us consider the MHD waves propagating inside the magnetic loops. The MHD waves consist of the forward wave Bfor 
and the reflected wave Bj-cc {B(ai and -Brof are the magnetic fields in the poloidal plane) . When the wave -Bjclk is transmitted 
to region III (see Fig. |8]), the electromagnetic energy can be converted to the kinetic and thermal energies in region III. The 
magnetic field Be given in equation (|6ip can be expressed by the superposition of the forward and reflected waves. When the 
density enhancement appears ahead of the magnetic loop in region II, the forward waves can be partially reflected by it. The 
condition for the perfect reflection should be determined by the wavelength A and the thickness of the density enhancement 
d ~ (a — b)t. This situation is analogous to the enhancement of the reflection rate by coating a glass with dielectric medium. 
The reflection rate becomes maximum when the width of the dielectric medium d satisfies d — {21 + l)A/4. When A — 27r/fc, 
this condition coincides with equation (|62p . Note that the parameter k in equation (|62p is not exactly the wave number but 
it determines the profile of the magnetic fields (see equations p9p . (|40p . and (|4ip for the definition of k). When the condition 
(|62|) is satisfied, the MHD waves propagating in the +r direction are totally refiected by the density enhancement produced by 
the loop expansion. For the dipolar solution, the magnetic energy is transmitted to r > R{t) because the density enhancement 
does not appear (see the right panels of Fig. |3]and Fig. |4}. 
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Figure 8. Schematic picture of the propagating waves. 
The forward wave -Bfor and the reflected waves -B^ef 
propagate inside the loops. Solid curves show the loop 
top at r = R{t), while the thin curves show the mag- 
netic field lines. The superposition of these waves de- 
termines Bg in region II. The leak wave Bioak appears 
ahead of the loop top in region III. 
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Figure 9. Contour plots of the flux function A/{Aqo?) of the shell solution for 
k = 7r/[2{a — b)\ (left panel) and k = 37r/[4(a — b)] (right panel). 

When k > 7r/[2(a — &)], the magnetic shell recedes from 77 ~ a to the region b < rj < a, and a flux rope appears 
around ri = a ahead of the magnetic shell. Fig. |9] shows the contour plots of the magnetic flux A of the shell solutions for 
k = n/[2{a — b)] (left panel) and k = 2,-k /[A{a — b)] (right panel). 



3.4 Flux Rope Solutions 

Flux rope solutions which include flux ropes inside the expanding magnetic loops are constructed by the flux function (|42 
By substituting equation (|42p into equation (|32p . the function Q can be written as 



sm 6, 



(1 



3rA-(7?) sin"e, 



(63) 



(64) 



where Qo.n is a constant and subscripts I and II denote region I and region II, respectively. The pressure and density functions 
(i.e., P and D) can be obtained by substituting equations (|42|) . (|63|l . (|64|) into equations (|31|l and (|33p . The functions P and 
D obtained in region I and region II are given in appendix [B] 

The pressure and the gas density consist of three parts, the isotropic part Po and parts representing the interaction 
with the electromagnetic force by the poloidal and toroidal magnetic fields. Pa and Pq (see equations (IBip and (|B7|) for the 
pressure and equations ()B4|I and (|B10|) for the gas density). 

Fig. 1101 shows the contour plots of the magnetic flux A (left), the poloidal part of the pressure Pa (centre), and that of 
the gas density Da (right) in rj/a — 9 plane. Fig. Illl shows contour plots of the toroidal magnetic fleld (left), the toroidal 
part of the pressure Pq (centre), and that of the gas density Dq (right) in -q/a — 6 plane. The parameters are taken to be 
a = 0.8, b = 0.75a, and m — n = 4 in both figures. The flux ropes exist behind the loop top (see the left panel of Fig. llOjl . 

The magnetic fields in region I and II are explicitly given by 



J 2Aoa^ 



B 



Bl = 



(65) 



(66) 
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rjkt vkt 

Figure 10. Contour plots of the magnetic flux A (left), the poloidal part of the pressure (centre), and that of the gas density Dj^ 
(right) for the flux rope solution in r]/a — 9 plane when a = 0.8, b = 0.75a. 



Br = ^ - ffis^ (68) 

bV = ^ ^ ^ 1 4^ [1 - ir/tf] sin^mrA))cos(T(./t)) _^ 1 

[l-{r/tf]i I ^ ^ ^ sin*(r(a)) W ;j , 

Here the subscripts / and II denote the magnetic fields in region I and region II, respectively. Similarly to the dipolar and 
shell solutions, the parameter m represents the Fourier modes which specify where the magnetic twist is injected. 

In the limit that t ^ r, the flux rope solution reduces to the dipolar solution given by equation (|52|) . (|53p and (|54p . In 
this limit, the magnetic field becomes stationary and radial. 

Note that the field component Br and Be are exactly zero but Be is not zero at r = R{t) unless the condition (|62p 
is satisfied. The electromagnetic energy is transmitted to r > R{t) unless the condition l|62p is satisfied as discussed in the 
previous subsection. When equation (|62p is satisfied, the Poynting flux is totally reflected at r = R(t) and the electromagnetic 
energy is not transmitted to r ^ R{t). 



4 PHYSICAL PROPERTIES 

Here we discuss physical properties of the three solutions we derived in § |3l In this section, we organize our discussion into 
four parts. First we consider the energetics. Second we show the shell and flux rope structures derived in I3.3l and 13.41 inside 
the magnetic loops. Third we study the relativistic effects, especially the role of the displacement current. Finally, we apply 
our solutions to SGR flares. 
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Figure 11. Contour plots of the toroidal magnetic field (left), the toroidal part of the pressure Pq (centre), and that of the gas 
density Dq (right) for the flux rope solution in rj/a — 6 plane when a = 0.8, 6 = 0.75a, and n = m = 4 . 



4.1 Energetics 

First let us consider the dipolar solution without the toroidal magnetic field (i.e., Qo,n = 0) for simplicity. Total energy £ 
contained inside the expanding magnetic loops is given as 



£ = K + Ui,, + Utb + UE + UM + W, 
where 



K = [ dVp-f\ 
Jv 



-L 



Ut^= / dV 



V 



r-1' 



Ue= I dV—, 



Um= dV 



8tt 



W = - dV 



(71) 

(72) 
(73) 
(74) 
(75) 
(76) 
(77) 



are kinetic, thermal inertial, thermal, electric, magnetic, and gravitational potential energies, respectively. Since the solutions 
we derived describe the freely expanding magnetic loops, i.e., Dv/Dt = 0, the total kinetic energy K given by 



P pa PZ7V PTT 2 7-) 

K= fyy^dV = dri dd d<l> ^ 
Jv Jo Jo Jo vl 



■ sm t>, 



(78) 



does not change with time. Here the total kinetic energy is integrated inside the spherical surface of r = R{t). Other energies 
can be evaluated by carrying out the integration directly. The non- kinetic part of the total energy £' = Uin + Uth+UE + Um+W 
contained inside r = R{t) is then given as 
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Toioida] magnetic field 



Figure 12. A schematic picture showing an expanding magnetic loop. Toroidal magnetic field is created inside the magnetic loop due 
to the twist injection from the surface of the magnetar. The magnetic pressure gradient force plus the pressure gradient force balances 
with the magnetic tension force by the poloidal magnetic field. 



Since the thermal, gravitational potential, and magnetic energies contain infinity due to the divergence of p, p, and Br 
at r = ( see equations (|45[) . (|49|) . and (|51|l . respectively), we renormalized the infinite parts of Uth, Um, and W to zero (see 



Low 198 



The non-kinetic part of the total energy £' depends on the amplitude of the poloidal magnetic field Ao, but is independent 
of the isotropic component (i.e., Po and Do). The isotropic component does not contribute to the total energy because the 
thermal energy of the isotropic plasma cancels with that of the gravitational potential energy. 

The energy £' diverges ai t = because we assumed a point mass at the origin. In magnetars, since the magnetar has a 
finite radius Rs, the self-similar expansion will take place when r > ro > Rs and t > to. Let us denote the total energy and 
the non kinetic part of the total energy inside the spherical surface of ro = R{to) as £q and £q, respectively. The expansion 
takes place when £o > £q. Since the total kinetic energy K = £o — £q does not change with time and £'{t) given in equation 
H79[) decreases with time for the dipolar solution, £{t) — £'(t) + K < £q + K = £q. The released energy £q — £{t) is carried 
away to r > R{t). This can be confirmed by integrating the energy conservation equation inside the spherical surface r = R{t) 



V(t) 



d_ 

dt 



, , 2 +B'^ GMpj 



8-K 



cPr+ I V- 



, , 2 ExB GMp-i 

(p + 4p)7 V H V 

A-K r 



d^r = 0. (80) 



'V(t) 

Since this integration is cumbersome, we do not show the details of the calculation. We have to point out that the integration 
with the volume V cannot be exchanged with the time derivative in the first term since the volume V changes with time. 

Next let us consider the case Qo,^ 7^ 0. Since the integration of equations (|72p - (|77p is complex, we evaluate the total 
energy inside the closed boundary by a different method. Note that equation p6ll . which represents equations of motion in 
the self-similar space, indicates that the self-similar equations we derived are closely related to the static relativistic MHD 
solutions except the existence of the thermal inertial term and the electric field. We can derive the virial theorem for the 
relativistic self-similar MHD (see appendix [Cj; 

3{r ^ l)Uti, + Uin + Um + Ue + W = H + S, (81) 
where 



J pr-dA-^ J {2[{r ■ E){E ■ dA) + {r ■ B){B ■ dA)] - {E'' + B^){r ■ dA)} 



(82) 



and 

5= Idvl^ 

dt 



Here A is the surface enclosing the volume V . The non-kinetic part of the total energy £' can be written from equation (|8ip 
as 
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Figure 13. k dependence of the non kinetic part of the total energy (solid curve) for the shell and flux rope solutions. 



= -(3r-4)l7th +-W + 5. (84) 

We can evaluate £' inside the expanding spherical surface oir = Rif) by using the fact that p = p = Br = = at r = _R(t) 
as 



3t ' 

16Aga^fc2(l-a^)cot^r(a) 

at ' 

WA^a^k^ cot"^ T(a) 



(dipolar solution), 

(shell solution), (85) 
, (flux rope solution), 

In all solutions, the non-kinetic part of the total energy does not depend on the toroidal magnetic field because the toroidal 
magnetic field does not change the dynamics of the expanding magnetic loops in the self-similar stage. This can be understood 
from the fact that when we take Qo.n = 0, the solutions we derived satisfy equations (|3ip . (|32|l and (|33p without any 
modification on the poloidal magnetic field. To understand this reason, let us consider the equation of motion in the 6 
direction. Since v = vBr, the force balance should be attained in the 6 direction (see Fig. I12|l . The pressure Pq is smaller for 
larger toroidal magnetic fields because Pq is proportional to — Qo,n (see equations (|46|l for the dipolar solution, (|A4[) and (|A11[) 
for the shell solution, and (|B3P and (|B9|l for the flux rope solution). For larger toroidal magnetic fields, the magnetic pressure 
gradient force by the toroidal magnetic field balances with the magnetic tension force from the poloidal magnetic field. As a 
result, the existence of the toroidal magnetic field modifies the plasma distribution, but does not change the dynamics. 

Fig. [13] shows the k dependence of £q for the shell and flux rope solutions (note that both solutions have the same k 
dependence). When the condition (|62|) is satisfied, the total energy contained inside the magnetic loops is equal to K and 
conserved for the shell and flux rope solutions because the energy flux is zero at r = R{t). When (|62|) is not satisfied, the total 
energy is larger than K by £' (see equation (|62p ). The excess energy is carried away to the region III (r > R{t)) to attain the 
free expansion, i.e., Dv/Dt = 0. 



4.2 Shell and Flux Rope Structures 



Let us examine the density and pressure distribution for the shell solution derived in § 13.31 We define the density and pressure 
enhancements as 

„4 

' Pa, (86) 



APa = 



2717?" 



Aga* sin^ 6 



ADa = 



APq = 



TiGMTfil-rf) 



2Ala* sin^ 6 



Da, 



A-K{m -I- n - 2)77^(1 - rf) 



nQo ,mQo,n sin" 



Pq, 



(87) 
(88) 
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Figure 14. The pressure enhancement AP^i (left) and the density enhancement ADa (right) for the shell solution are shown as a 
function of ri/a. Solid curve is for 6 = 9a/10, while the dashed and dot-dashed ones are for b = 3a/4 and 6 = a/2, respectively. Other 
parameters are fixed at a = 0.8 and k = 7r/[4(a — b)\. 



ADo = 



2nGM{m + n- 2)r?(l - rj'^f 



nQo,mQo Sin 



m-\-n — 2 . 



Dq, (89) 



These functions are normalized to be unity in region I, where the poloidal magnetic field lines are radial for the shell 
solution. In Fig. 1141 the pressure and density enhancements, APa and ADa are plotted for b = a/2, 3a/4, 9a/10 when 
a = 0.8 and k = 7r/[4(a — b)]. In all three cases, the pressure and density pulses appear at the top of the magnetic loops. 
Their amplitudes are larger for a thinner shell. The peak of the pressure enhancement appears behind that of the density 
enhancement. This structure comes from the requirement for the force balance with the gravity. As mentioned in 14.11 this 
relativistic self-similar solution is similar to the static solution in which the force balance is attained. As plasma is swept up 
into the shell, the density increases inside the shell. To support the gravity by this excess density, the pressure gradient appears 
behind the density enhancement. The density decrease behind the pressure enhancement also comes from the requirement for 
the force balance. Since the decrease of the density enables the buoyancy force to push the plasma in the radial d irection, th is 



buoyancy force maintains the pressure pulse. These structures are identical to those in non-relativistic solution (|Lowlll982l ). 

Fig. US] plots APq and ADq for b = a/2, 3a/4, 9a/10 when a = 0.8 and k = n/[4{a - b)]. As mentioned in § |31 the 
Lorentz force exerted by the toroidal magnetic fields always reduces the pressure. A local minimum of the density enhancement 
ADq locates behind a local maximum of dAPg /drj. This structure also comes from the force balance. Pressure gradient force 
balances with the buoyancy force in the rarefied region. 

Next we examine the structure of the flux rope solution derived in § 13.41 We define the normalized toroidal magnetic field 
strength as 

^B4V) = ^ ^"'Lr, - (90) 
Qo,n sm" 6 

Solid curve in Fig. 1 16 1 shows AB<^ as a function of r; for b = 0.95a, while the dash and dot-dashed ones show that for b = 0.8a 
and b = 0.65a, respectively. Other parameters are fixed at a = 0.8 and k — 7r/[4(a — b)]. The toroidal magnetic field has a 
peak inside the fiux rope. Its amplitude is larger for a larger a and a thinner shell. The shell structure also appears behind 
the loop top (see Fig. 1101 and [TT|) . Solid curve in Fig. 1171 shows AB^, which corresponds to the magnetic pressure by the 
toroidal magnetic field, as a function of rj for a = 0.8, b = 0.95a, and k — 7r/[4(a — b)]. Dashed and dot-dashed curves show 
APq and ADq, respectively. Plasma density decreases inside the shell. The decrement of the plasma density leads to the 
buoyancy force which balances with the pressure gradient force in front of the shell. Behind the shell, the pressure gradient 
force balances with that of the magnetic pressure. This effect is more prominent for the flux rope solution than for the shell 
solution since the magnetic pressure is enhanced inside the fiux rope. 



4.3 The Role of the Displacement Current 

We showed that Pq is always negative. On the other hand. Pa can have either positive or negative values. In this subsection, 
we obtain the condition for Pa < 0. 

First let us consider the dipolar solution derived in § 13.21 The condition that Pa given by equation (|45p is positive in 
^ ?7 ^ a is given by 
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Figure 15. The pressure decrement APq (left) and the density decrement ADq (right) by the azimuthal field are shown as a function 
of Tj/a for the shell solution. Solid curve is for b = 9a/10, while dashed and dot-dashed curves are for b = 3a/4 and b = a/2, respectively. 
Other parameters are fixed at a = 0.8 and k = 7r/[4(a — b)]. 



2a - -ia^ - T?"* + Irf" ^ 0. (91) 
The critical value a* for Pa >0in0S7? = ais 



a. = V69+p^ , 0.958. (92) 

When a > a*, Pa has negative values in the domain ^ 77 ^ a. Since a denotes the expansion speed of the magnetic loops at 
r = R{t\ the above condition indicates that Pa can be negative for faster expansion. 
Next let us calculate the azimuthal component of the current density, 

U = >ot + jdisp, (93) 
where 

Jrot = = - — IT 5 sme, (94) 

1 dE^ Aoa^ rf(G-3a^-5rf + 2rf) . 

jiisp = -- T7r = --, — TT ^ 5 ^sm6l, (95) 

An dt Attt^ (X_jy2)§ 

and 

^'-1^ JT^i ^''^ 

The current jrot is always positive, while jdisp has negative values for a larger a in ^ 77 ^ a. The displacement current jdisp 
cannot be ignored for a larger a and it reduces the azimuthal current j^. Thus the current changes its sign for a larger a. 
Remember that the pressure P is determined by the d component of the equation of motion given by 

(-Vp + j X B + p,B), = 0. (97) 
According to the definition of Pa and Pq, the poloidal component of equation (|97|l is given by 

i^-i.^^- (98) 

Since both dPA/dO and j0-Br depend on 9 by sin S cos 6, and Br/ cos is positive, the sign of Pa is determined by that of 
j^/sinS. Thus Pa can be negative when the displacement current jdisp dominates the current j^ot- The condition that j^ ^ 
coincides with the condition that Pa ^ (i.e. a ^ a*, where a* is given by equation (HU). 

Next let us consider the shell and flux rope solutions. Since these solutions are more complex, the equation Pa = is 
solved numerically. Instead of using parameters a and b, we introduce the following parameters, 

Knax = a, (99) 
A=^. (100) 
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Figure 16. Distribution of AB^ for the flux rope solution as a 

function of ri/a for a = 0.8. Solid curve is for b/a = 0.95, while Figure 17. Distribution of AB^ for the flux rope solution as a 

dashed and dot-dashed ones are for bja = 0.8 and b/a = 0.65, function of -q/a for a = 0.8 and b/a = 0.95 (solid curve). Dashed 

respectively. and dot-dashed ones denote APq and ADq, respectively. 

These parameters denote the maximum speed of the expanding loops and thickness of the shell behind the loop top, respec- 
tively. Fig. [18] is a diagram showing whether the solution that Pa ~ exists for the shell and flux rope solutions in the 
parameter space of Vinax and A for k = n/[4{a — b)]. Solutions where Pa = exist in the shaded area for the shell solution 
and in the grey area for the flux rope solution. Similarly to the dipolar solution, the effect of the displacement current is more 
prominent for a larger Knax. Generally, Pa is smaller for a larger Vinax and thicker shells in parallel that the displacement 
current becomes important for larger Vinax and A. For the flux rope solution, the displacement current is important not only 
in region II but in region I (right bottom region in Fig. [18}. In this case. Pa is negative for b ^ 0.817. 



4.4 Application to SGR Explosions 



SGR flares can be triggered by energy injection into magnetic loops at the surface of a strongly magnetized neutron star (e.g., 
[Lvutikov,. 20061 ) . When sufficiently large energy is injected, the magnetic loops will become dynamically unstable, and expand 
relativistically. Magnetic energy release in the expanding magnetic loops can be the origin of SGR flares. The expanding 
magnetic loops will also produce magnetosonic waves propagating ahead of the loops. High energy particles can be produced 
in the magnetic reconnection inside the loops, and in shock fronts formed ahead of the loops. 

In this paper, we did not solve the structure of the region ahead of the magnetic loops (r > R{t)). When the outer 
region is a vacuum, electromagnetic wave s will be emitted from th e boundary at r = R{t). When the plasma density is much 
larger than the Goldreich- Julian density ( Goldreich fc Julian 19691 ) and the wave frequency is much smaller than the plasma 
frequency, the outer plasma can be studied by using MHD equations. It will be our future work to connect the self-similar 
solutions inside r = R{t) and the solutions in r ^ R(t). 

Now let us es timate the energy for th e SGR explosion based on the self-similar solutions. Let us take the field strength 
to be 10^^ Gauss jKouveliotou et al.|[l998l : llbrahim et al.ll20o3 : llbrahim et al.ll2003h at the stellar radius Rs = lO'' cm. This 
leads to 



'^^'^^ — iQ^^ Gauss. 
Ri 



(101) 



By assuming that the self-similar expansion begins when to ~ Rs/a, the released energy from the expanding magnetic loops 
can be estimated from equation (|85p as 



(dipolar solution), 



6 X lO"*® erg, 
2 X lO'^^'A"^ erg, (sheU solution), 
8 X lO^'^A"^ erg, (flux rope solution). 



(102) 



Here we take a = Vmax — 0.9 and fc = 7r/[4(a — 6)1. Th ese results agree with the observed energy of SGR giant flares 
( Hurlev et al. 2005 : Palmer et aLlboOsI : Terasawa et al.ll2005 ). Note that the total energy contained in the expanding magnetic 
loops is more energetic for thinner shells. The non-kinetic part of the total energy £' is inversely proportional to the square 
of the shell thickness. When some fraction of the kinetic energy K is converted to the electromagnetic energy, the released 
energy can be larger than that estimated by equation (|102p . 
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Figure 18. The diagram showing where the solution = exists in the parameter space Vniax and A = (a — b)/a when k = 7r/[4(a — 6)] 
for the shell solution (the shaded area) and for the flux the rope solution (grey area). 



5 SUMMARY & DISCUSSIONS 

By extending the self-similar solutions derived by Low ( 19821 ). we derived self-similar solutions of relativistically expanding 



magnetic loops taking into account the toroidal magnetic fields. The dipolar solution derived in 13.21 gives us an insight into 
the relativistic expansion of the magnetic loops because of its simplicity. However, the shell and flux rope solutions derived 
in § 13.31 13.41 have more physically interesting properties such as an enhanced magnetic pressure at the shells and flux rope 
structures. Such configurations might be more probable for SGR flares. 

The equations of motion in the self-similar stage are similar to those of the static equilibrium state except the existence 
of the relativistic thermal inertial term and the electric field. This fact allows us to evaluate the non-kinetic part of the total 
energy in the magnetic loops by using the virial theorem. The magnetic loops with shell or flux rope structures carry more 
energy than the simple dipole solution. The energy is comparable to the observed energy of the SGR giant flares. 

In relativistically expanding magnetic loops, the effect of the displacement current becomes important. In dipolar solution, 
the displacement current becomes larger than the real current V x i?/(47r) for faster expansion speed (Knax > a*)- This effect 
reduces the toroidal current and weakens the magnetic tension force. To balance the reduced magnetic tension force, the 
pressure decreases. 

We found that the energy is transferred to r > R{t) in dipolar solutions. In the shell and flux rope solutions, the energy 
is transferred to r > R{t) unless condition (|62|l is satisfled. The condition can be interpreted as that for the total reflection 
of the MHD waves in the shell. Dipolar solutions always have leakage (transmission of Poynting flux to the region r ^ R{t)) 
because Bg ^ at r = R{t). The shell and flux rope solutions have perfectly reflecting solutions in which the total energy in 
r < R(t) is conserved. It means that the solutions are energy eigenstates of the system. The eigenstates can be obtained by 
adjusting the parameter k. 

In this paper, we obtained solutions for freely expanding magnetic loops, i.e., Dv/Dt — 0. We assumed that the magnetic 
loops have sufficiently large energy to drive the expansion. When the flux function A increases with time, the toroidal magnetic 
fields will also increase with time. The toroidal magnetic fields will then affect the dynamics through the magnetic pressure. 
Such solutions can describe the accelerating magnetic loops. 

Magnetic fields can be expressed as the sum of the Fourier modes in the polar angle. The modes and their amplitudes 
should be determined at the boundary where the magnetic twist is injected on the surface of the star. It is not shown but 
we can construct more complex solutions that the poloidal magnetic fields are expressed by the sum of the Fourier modes, 
i.e., A (X sin" 6. In actual explosion, the opening angle of the expanding magnetic loops depends on the location at which 
the magnetic twist is injected on the surface of the central star. Such a solution may be expressed as the sum of the Fourier 
modes for the poloidal and toroidal magnetic fields. We should note that SGR fiares are not necessarily axisymmetric. Models 
including the non-axisymmetrically expanding magnetic loops will be a subject of future works. 
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APPENDIX A: CONSTRUCTION OF THE SHELL SOLUTIONS 

The functions Q and P in region I are given by 

n 

P'{v,e) = P„(r?) +Pi(r,,e) +P<^(r;,e), (A2) 
where P4 and Pq are given by 
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(A3) 



Pq{v,0) = < 



E 



nQo,mQo,7i 



m + n — 2 / 



47r(m + n — 2)772(1 — rj^) 



nQo 



log (sine), 



for m + n 7^ 2, 
for m + n = 2. 



m + n— 2 

The functfon D in region I is given by 
D\r^,e) = i3o(r?) +i3i(r?,e) + 0^(r;,e), 
where and Dq are given by 



7rGMr?3(l_^2) 



I>q('?,^) 
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m + 715^2 
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2-KGMr){l - 772) 



■log(sin6i), 



for 771 + 7T, 7^ 2, 

for m + 71 = 2. 



m + n — 2 

The functions Q and P in region II are given by 
Q"(r,,e) = ^^At(7,)sin"0, 

n 

where and Pq are given by 



Pi'iv, 
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for 7n + n 7^ 2, 

for 771 + 71 = 2. 



The function D^^ is given by 

D"(v,e) = Doiv) + d'J{v,o) + Dlf{n,e), 

where D^^ and Dq are given by 
2y 
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for 771 + 71 = 2. 

Here 2(77) is a function of 77 given by 

=(77) = l-fc77(l + 7?^)(l- 377^)cotr(7?)+fc^7?^(l-77^)(l + 377^) [1 - 3cot^r(77)]+fe^77^(l-77^)^ [3cot^r(r7) - 5cotr(7?)] 
Poiv) is an arbitrary function of 77 and is related to the function Doirf) through equation (|48|) . 
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The function P in region I is given by 
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where and Pq are 
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The function D in region I is given by 
Z)^(r?, 9) = Doiv) + Diirj, 9) + D'q{71, 9), 
where and Dq are 
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The function P in region II is described as 
P"(7?, 0) = Po(r?) + Pr(7?, 9) + P^'in, 9), 



for 7Tl + 71 2, 
for TTl + 71 = 2. 



where P" and Pq are 
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The function D is given by 

D" (77, 9) = Do (77) + D'J (77, e) + D'q' (v, ( 



where D'J and Dq are 
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The functions do, di, d2, da, and Y are given by 
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The functions Pq and Do are related each other through equation (|48 



APPENDIX C: THE VIRIAL THEOREM OF THE SELF-SIMILAR RELATIVISTIC MHD 

The virial theorem in non-relativistic MHD was derive d by Chandrasekhar fc Fermi ( 1953h . Low (1982) apphed it to the 
expanding magnetic loops by evaluating the surface term. Landau fc Lifshit j ( 1975 ) derived the theorem for a relativistic case 
in elegant way by integrating the energy momentum tensor. In this appendix, we derive the virial theorem for a relativistic 
self-similar MHD. 

We start from the equations of motion in self-similar stage given by (|36[) . Taking the inner product with r and integrating 
it within a volume V, we obtain 

F p7^ii^(r • er) 



dV 



r -Vp + r ■ {peE +j X B) 



^(r-.e.) 



: 0. 



(CI) 



The first term is the thermal inertial term Uh-i and the forth term is the gravitational potential energy W. Integrating the 
second term by parts, we obtain 



dVr ■ Vp = -3(r - l)(/th + pr dA 



(C2) 



where ^ is a closed surface of the volume V. The third term can be rewritten by using the Maxwell equations as follows, 

d 



dV[r-[p,E+3xB)] = - / dV-[r-S) 



dVr ■ (V ■ 



(C3) 



where S and cr are the Poynting flux and the Maxwell's stress tensor, respectively. Note that the time derivative cannot be 
exchanged with the integration with volume in the first term since the volume V changes with time. 
The second term on the right hand side of equation (|C3|I has a form 

(C4) 



j dVr ■ {V ■ <t) ^ + Um + ^ j {2 [(r ■ E){E ■ dA) + (r • B){B ■ dA)] - {E^ + B^){r ■ dA)} 



where Ue and Um are the electric and magnetic energies given by equations (|75p and H76p . respectively. By using these results, 
we obtain the virial theorem for the relativistic self-similar MHD; 



3(r - i)c/th + Uin + UM + UE + w = n + s, 

where 



H = J pr-dA - ^ J {2[{r ■ E){E ■ dA) + {r ■ B){B ■ dA)] - (E'' + B^){r ■ dA)} , 



and 
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4-K 



(C5) 
(C6) 

(C7) 



Here, K is the kinetic energy given by equation 

Readers may wonder why the thermal inertial term appears. Actually, the kinetic, thermal, and thermal inertial energies 
should not be considered separately because they depend on the frame of reference. Even so, we used this definition through 
the paper to make clear the difference between the non-relativistic and relativistic expansions. We can easily figure out that 
the total plasma energy can be described as the sum of these energies, as 



K + f/th + Uin 
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p 7 
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